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EXISTENCE AND STABILITY OF SOLITONS FOR 

FULLY DISCRETE APPROXIMATIONS OF THE 
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^^ Abstract. — In this paper we study the long time behavior of a discrete ap- 

(-y-j proximation in time and space of the cubic nonhnear Schrodinger equation on 

^Nj the real Une. More precisely, we consider a symplectic time splitting integrator 

applied to a discrete nonlinear Schrodinger equation with additional Dirichlet 

I I boundary conditions on a large interval. We give conditions ensuring the ex- 

■^^ istence of a numerical soliton which is close in energy norm to the continuous 

^7' soliton. Such result is valid under a CFL condition of the form rft^^ < C 

• where r and h denote the time and space step size respectively. Furthermore 

T^H we prove that if the initial datum is symmetric and close to the continuous 

5^ soliton 77 then the associated numerical solution remains close to the orbit of 

C^ r;, r = Ua{e*"r7}, for very long times. 



> 

cn 

Q^ 1. Introduction 

^^ We study numerical approximations of solitons of the focusing nonlinear 

CO Schrodinger equation (NLS) on the real line: 

(N (1-1) ^V'i = -i^xx - IV'I V, a; G M, t G M. 

• • This equation is a Hamiltonian partial differential equation (PDE) associated 
with the Hamiltonian function 



;l2) H(:4,) :-- 



w.P-^ 



dx, 
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and preserving the L^ norm 

(1.3) iV(^) := / IV'l^dx. 



The goal of this paper is to understand the long time behavior of numeri- 
cal integration algorithms for initial data close to the solitary wave solution 
V'(t,x) = e*'^*r/(x) where 

1 f X 

(1.4) r]{x) := -/|sech y- 

and A € M is the Lagrange multiplier associated with the minimization of H 
under the constraint A^ = 1. It is well known, see for instance [20, 13, 14, 12] 
that this solution is orbitaly stable in the sense that for a small pertuba- 
tion of the initial data, the exact solution remains close to the orbit of t] for 
all times. Here we will only consider symmetric initial conditions satisfying 
ip(x) = 'tp{—x), a property that is preserved by the flow of (1.1). In this set- 
ting, the orbital stability of the continuous soliton can be described as follows: 
Let 

(1.5) T:=[j{e^-r^(x)} 

and assume that ip{0, •) is a symmetric function satisfying dist('0(O, •),r) < 6 
for some 5 sufficiently small, then for all times t > 0, if ip{t, •) denotes the 
solution of (1.1), we have 

(1.6) Vt>0, dist{%b{t,-),r) <C6, 

where C is a constant independent of 5 and t, and where the distance is mea- 
sure in H^ norm. The present paper deals with the persistence of this result 
by fully discrete numerical methods. It is an old problem that was pointed out 
in several papers in the last 30 years, see for instance [7, 18, 8, 5], and the 
numerical approximation of (1.4) over long times has now become a classical 
benchmark to test the performance and stability of numerical schemes, see for 
instance [1, 9, 4] and the references therein. However, as far as we know, 
no result of the form (1.6) has been proven in the literature for fully discrete 
approximations of (1.1) (see however [2, 5] for the space discretized case). 

In particular, the effect of the time discretization yields many mathematical 
difficulties. Duran & Sanz-Serna gave in [8] some asymptotic expansion of 
the numerical solution close to a soliton, but the lack of a modified energy 
acting on H^ and preserved over long time by the numerical scheme (the so 
called backward error analysis) was an obstruction to define a possibly stable 
numerical soliton. Here, we take advantage of a recent construction of such a 
modified energy given by Faou & Grebert in [11] to show the existence and 
stability of a modified soliton that is close to (1.4) in energy norm. 
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In this paper, the discretization of (1.1) we consider are made of three levels 
of approximations: 

— A space discretization, where we use a grid with mesh size h > 
made of an infinite collection of equidistant points of M. The equation 
(1.1) is then approximated by the discrete nonlinear Schrodinger equation 
(DNLS) where the Laplace operator is replaced by its finite difference 
approximation over three points. 

— A Dirichlet cut-off, where we replace the integrability condition at 
infinity of the derivative of "0 by a Dirichlet boundary condition at the 
boundary of a large window of size 2Kh where K >> 1. 

— A time discretization algorithm to integrate the DNLS equation 
with Dirichlet boundary condition. This discretization introduces a last 
parameter r which represents the time step. To do this we consider a 
symplectic time splitting algorithm where the kinetic part and potential 
part are solved alternatively as described for instance in [19]. 

Each of these three levels of discretization relies on discretization parameters. 
In this paper, we prove orbital stability in the sense of (1.6) for the numerical 
solution, where the distance to T is estimated in terms of the three discretiza- 
tion parameters h, K and r. 

We first present some numerical experiments showing that the solitary wave 
rapidly disappears if either the algorithm of integration is not symplectic, or 
if it is symplectic, but used with a too large CFL number r/i~^. 

The proof is organized as follows: we first recall in Section 4 the main 
arguments of the proof of the orbital stability result in the continuous case, 
following in essence the presentation made in [12]. We then give in Section 5 
an abstract result showing that if the energy space H^ is well approximated 
by the space discretization, and if the numerical scheme preserves - or almost 
preserves - modified L^ norm and energy functions that are close to the exact 
ones, we can obtain orbital stability results with precise bounds depending on 
the parameters. We then apply this formalism in Section 6 to our three levels 
of discretization. 

As the proof of orbital stability result is based on the variational charac- 
terization of the solitary wave and thus heavily relies on the preservation of 
the energy and L? norm, long time bounds can be straightforwardly obtained 
for energy and L? norm preserving schemes such as the Dufour-Fortin-Payre 
scheme, see [7]. This follows directly from the analysis of the space discretized 
case (see also Remark 2.6). 

The cornerstone of the analysis of splitting method is the construction of 
the modified energy. Recall that in the finite dimensional case, the existence 
of modified energy is guaranteed by Hamiltonian interpolation: see [3, 15, 16] 



4 DARIO BAMBUSI, ERWAN FAOU & BENOIT GREBERT 

but cannot be applied straightforwardly to Hamiltonian PDEs unless unrea- 
sonable a priori assumptions are made on the regularity of the numerical so- 
lution, which prevents a fair use of the bootstrap argument underlying the 
orbital stability methodology. Here we take advantage of the recent backward 
error analysis result of [11] to construct a modified energy acting on H^ for 
splitting methods applied to (1.1). Actually we give a simplified proof of a 
simpler version of the result presented in [11] or [10], which has some interest 
in itself. 

Using this result, we then prove an orbital stability result for fully discrete 
splitting method applied to (1.1) with a CFL restriction, and over very long 
times of the form nr ~ t~^ , where M > is an integer number depending 
on the CFL. 

2. Three discretization levels and main results 

We now describe more precisely the three levels of approximation of (1.1) 
mentioned in the introduction. At each step, we state the orbital result that 
we obtain. 

2.1. Space discretization. — Having fixed a positive parameter h we 
discretize space by substituting the sequence tp£ ~ tpihi), £ G Z for the function 
ip{x), and the second order operator of finite difference A^ defined by 

(2.1) [Ahi^je ■■= p , 

for the Laplace operator —dxx- The NLS is thus reduced to the discrete 
nonlinear Schrodinger equation (DNLS): 

(2.2) ii;e = -^(V-^+i + ^i-i - l^i) - [V-^l V^ ^ e ^ . 

where t i— )• ^(t) = (^^(t))^^^ is an application from M to C^. With this 
equation is associated a Hamiltonian function and a discrete 1? norm given 
by 

2 I ,, |4' 



(2.3) Hu{^\^) = hY, 



^3 - V'i-1 



h 



IV', 



and Nh{i)) = hY,\^,?- 



The discrete space of functions is 
equipped with the discrete norm 
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Following [2], we identify Vh with a finite element subspace of H^{M.; C). More 
precisely, defining the function s : M — t- M by 



(2.4) six) 



if |x| > 1, 

x + l if -1<2;<0, 

-x + l if < X < 1, 

the identification is done through the map ih '■ Vh ^ H^CR; C) defined by 

(2.5) {^,}^.g^ ^ {ihi^){x) ■.= Y,^j [I - j) . 

jez 

Recall that T is the curve of minima of the continuous Hamiltonian and is 
given by (1.5). With these notations, we have the following result 

Theorem 2.1. — There exist 5q and ho such that for all 6 < 5o and h < ho, 

if {ip^)j(=,z £ Vh is such that 

dist(4V°,r)<(^, 

where the distance is measured in the continuous H^(K;C) norm, then the 
solution {tpj{t))j,zz of (2.2) satisfies 

Vt>0, dist{ihi'{t),r)<C{5 + h) 

for some constant C independent of h and 5. 

Notice that the DNLS flow is not defined globally everywhere, i.e. for all 
initial data in Vh and all times t. However since a solution of DNLS issued from 
an initial datum close to T remains unconditionally close to F, such solution 
is automatically global. 

2.2. Dirichlet cut-off. — In order to come down to a finite dimensional 
system we fix a large number K > 1, substitute the sequence —K,...,K for 

the set Z in (2.2), and add Dirichlet boundary conditions ■0-A'-i = '4'K+i = 0. 
The equation we consider is thus the (large) ordinary differential system 

f iipe = 2('0m + V'£-i-2V'^)-|V'dV^, -K<i<K 

(2.6) < ^J'^ 

( 4'±{K+i) = 0. 

Note that here, we use the convention that ^l^e = for all \i\ > K + 1, so 
that the previous system is indeed a closed set of differential equations. The 
corresponding discrete function space is 

(2.7) Vh,K-={ii^j)jez&Vh\i;j = for \j\>K + l}, 

on which we can define the Hamiltonian function and discrete L^ norm Hh^x '■= 
HhlvhK ^^'^ ^h,K '■= ^h\Vh K ^^ restrictions of the functions (2.3) to Vh^K C 
Vh- Similarly, we define ih,K = '^hWuK- ^"^ ^^^ following, we often use the 
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notation {'4'j)j=-k t° denote an element of Vh^x with the imphcit extension 
by for \j\ > K + 1 to define an element of (2.7). With these notations, we 
have the following result: 

Theorem 2.2. — There exist constants Ci, C2, Sq and cq such that for all 
6 < 60 and all h and K such that h + j^e~'^^^^ < eo, if {'ip'j)f=-K ^ ^h,K is 
such that 

dist(4,i^V°,r)<5, 

then the solution {il^j{t))j^_j^ of (2.6) satisfies 

Vt > 0, distiih,Ki^it),T) < C2 (5 + /i + J^e-^'^^ 

Remark 2.3. — The exponentially small term in the previous estimate rep- 
resents the effect of the Dirichlet cut-off. As we will see below, it directly 
comes from the fact that the function r/ is exponentially decreasing at infinity. 

2.3. Time discretization. — In this work the time discretization of (1.1) 
that we consider is a splitting scheme: we construct "0" the approximation of 
the solution il^{t) of (1.1) at time nr iteratively by the formula 

where the flow <I>p is by definition the exact solution of 

#^ = -IV'^PV^ (. = -K,...,K , 

in Vh^K which is given explicitly by formula <l>p('0)^ = exp(zT|?/'^p)'0£. The 
flow $^, is by definition the solution of 

(2.8) iV'£ = -p(V'm + V'^-i-2V'£), t = -K,...,K, 

with the convention -0^ = for |^| > K + \. The implementation of this 
numerical scheme requires the computation of an exponential of a tridiagonal 
matrix at each step. It could also be done in discrete Fourier space in which the 
operator on right-hand side is diagonal. The main advantage of this splitting 
method is that it is an explicit and symplectic scheme. 
Our main result is the following 

Theorem 2.4- — There exist constants Ci, C2, <^o o.'^d eq such that for all 
5 < 5q and all h, r and K such that h + -^e'^^^^ < eo and the following CFL 
condition is satisfied 

(2.9) (2M + 3)p<y, 
then if {'ip^)f^_i^ € Vh,K is such that 

dist(z;,,^V°,r)<(5, 
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we have 

(2.10) Vnr < t'^'^ , dist(4,^($^o$- ) V,r) < C2{6+h+'j^ + ^e-^'^'' 

Remark 2.5. — In the last estimate (2.10), the term r/h represents the error 
induced by the modified energy constructed with the method of [11] (see 
Section 7 below). Note that under the condition (2.9), this term is actually of 
order 0{h). 

Remark 2.6. — An alternative time approximation of (2.6) is the modified 
Crank-Nicolson scheme given by Delfour-Fortin-Payre see [7, 17] defined as 
the application ■0"' i— )• ■0""'"-'^ such that 

v^r' = V'? + y (A,(V'"+i + r))^ + jdv^r'i' + i^rT)(V'r' + ^d, 

for (. = —K, . . . ,K. It can be shown using a fixed point argument that for r 
sufficiently small, ^"+^ is well defined, and that this scheme preserves exactly 
the discrete L^ norm and discrete energy (2.3). Using this property, it can eas- 
ily be shown that the conclusions of Theorem 2.2 extends straightforwardly to 
this specific fully discrete case. Notice that this method has the disadvantage 
to be strongly implicit. 



3. Numerical experiments 

In this section, we would like to illustrate the results given in Theorem 2.4, 
and prove that if the CFL condition (2.9) is not satisfied, the stability estimate 
(2.10) is no longer true. In contrast, we show that if the CFL number is small 
enough, a numerical stability can be indeed observed. On the other hand, we 
show that for non symplectic integrators, even used with a very small CFL 
number, numerical instabilities appear. 

In a first example, we take h = 0.1875, if = 80 (so that Kh = 15), r = 0.2 
and the initial condition (1.4). The CFL number is equal to 5.7. We consider 
the integrator <I>^ o <I>p defined above. As mentioned in the previous section, 
the flow of ^pT can be calculated explicitely, while the computation of ^a - 
see (2.8) - is performed using the expm MATLAB procedure. 

In Figure 1, we plot the absolute value of the fully discrete numerical solu- 
tion ip"^ = ($^ o $p)"'('0'^). We can observe that the shape of the soliton is 
destroyed between the times t = 100 and 200. 

In a second example, we take the same initial data and parameters if = 80 
and h = 0.1875, except that we take a much smaller r = 0.001 making the 
CFL number equal to 0.028. However, we break artificially the symplecticity 
of the integrator by replacing the exact evaluation of the exponential in the 
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CFL = 5.7 




_ 0.5 WW s °-^ /'^^ 



Figure 1. Instability for t//i^ = 5.7 

flow ^^ by its Taylor approximation of order 2: 

exp(r^) ~ / + rA H ^^. 

As before, we observe in Figure 2 some instability phenomenon after some 
time, despite the fact that the CFL number is very small. Such an instability 
is due to the non symplectic nature of the integrator, which prevents the 
existence of a modified energy preserved by the numerical scheme. 

Finally, we consider the same initial condition and numbers K and h, but 
we take r = 0.02 making the CFL number be equal to 0.57 and we compute 
the exponential exactly making the scheme symplectic. In Figure (3) we can 
observe that the soliton is preserved for a very long time, up to t = 10® which 
corresponds to 2.10^ iterations. This result illustrates our Theorem 2.4. 



4. The continuous case 

Before giving the proofs of the Theorems presented above, we recall here 
the main lines of the proof of the orbital stability result in the continuous and 
symmetric case obtained first by [20] (see also [13, 14, 12]). The proofs of the 
discrete results will be essentially variations on the same theme. The method 



CFL = 0.028 



-15 -10 -5 
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_ 0.5 




Figure 2. Instability of non symplectic integrators 



is based on the variational characterization of the soliton r] as the unique real 
symmetric minimizer of the problem 



(4.1) 



min H{ih) 
Af(V)=l 



where H is the Hamiltonian (1.2) and N the norm (1.3). 

Remark 4-1- — By the method of Lagrange multipliers there exists A > 
such that 



-dxxr] - V 



-At/. 



Remark 4 •2. — We only consider the case where N{r}) = 1 in order to avoid 
the introduction of a supplementary parameter. It is clear to the reader that we 
could also consider the numerical approximation of any given soliton, provided 
that its L^ norm enters into all the constants appearing in the estimates below. 



In the following, we set 
We also define the real scalar product 



ijj{—x) = ip{x) }. 



(93,^)= Re / ip{x)tlj{x)dx. 



10 
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CFL = 0.57 
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Figure 3. Long time stability for t/K^ = 0.57 



This scalar product allows to identify //^(M;C) with the product /7-'^(M;]R) x 
H^{R; M) as follows: liip = -j={q + ip) and (p = -y^{q' + ip') where p, q, p' and 

q' are real symmetric //^(M;M) functions, then we have 



{'P,i^) 



{x)q'{x) + p{x)p'{x)dx. 



The real scalar product on H'^{R;C) ~ H'^{R;R) x H'^{R;R) is then given by 

(99,-0) = ((/7,'0) + {d^(p,dxi)), 
and we set 

for Lf = '^—T^- In the rest of this paper, we often amalgamate the two complex 
and real notations. 

In the following, we set 

(4.2) U{R) = yeV\ dist(^, T) < R}, 

where T is defined in (1.5), and the distance is measured in H^ norm. 
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Note that the Hamiltonian function H and the norm A^ are smooth in 
H^ (using the fact that H^ is an algebra). Moreover, these functions are 
gauge invariant, in the sense that for ah ip G H^ and all a S M, we have 
H{e'^"(p) = H{ip) and N{e^°'ip) = N{ip). Due to this invariance, it is immediate 
to realize that the whole manifold T is formed by minima of the minimization 
problem (4.1). Then it is well known [20, 13, 14, 12] that these minima are 
nondegenenerate in the directions transversal to the orbit F defined in (1.5), 
for symmetric functions. 

More precisely, following [12], we define the following set of coordinates in 
the vicinity of F: set 

(4.3) W = {u£V\{u,7]) = {u, iri) = 0}, 

equipped with the H"^ norm induced by the space V . As ir] is tangent to 
the curve F and orthogonal^^^ to ry, the previous W can be interpreted as the 
space orthogonal to the plane containing the planar curve F. Note that W 
is invariant under the multiplication by complex number: for any z G C, if 
u ^W then zu £ W. 

We define the map x as follows: 

(4.4) TxRxW 3 {a,r,u)^ x(a, r, u) = e*"((l + r)r] + u) e V, 

where T = ]R/(27rZ) is the one-dimensional torus. 

The following Lemma can be found in [12, Section 5, Proposition 1]. In 
our symmetric situation, we give here an independent proof that will later be 
easily transfered to the situation of discrete systems: 

Lemma 4- 3. — There exist constants rg and R such that the application x 
is smooth and bounded with bounded derivatives from T x [— ro,ro] x B{R) 
to V, and such for all ip G IA{K), there exists {a,r,u) G T x M x VF such 
that if = x(«)^;W)- Moreover, the application x~^ is smooth with bounded 
derivatives onlA{R), and there exists a constant C such that for all ip £ hl{R), 
we have 

(4.5) ||n(V)||^i<C7dist(V',F). 

Proof. — The first part of this lemma is clear using the explicit formula for 
X- To prove the second one, let us consider the projection of ifj onto the plane 
generated by {r],ir]): 

{ip, r])r] + {ip, ir])ir] =: z{ip)r] 

with z{^p) = {ip, rf) + i^ip, ir]) = J Tpfj £ C Note that the application tp i— )- z{ip) 
is smooth with bounded derivatives from V to C. Moreover, we have 

dist(V',F)2 > inf A^(V' - e^°??) > ||z(V')P - 1|. 



"^^^Recall that here ( • , • ) is a real scalar product. 
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Hence for R < 1/2 and for all il) G U{R), we have \z{'il^)\ G [1/2,3/2]. This 
shows that the apphcations 

U{R) 9 V ^ a(V') = arg(z(i/')) G T 
and 

U{R) 9 V ^ ^W = kWI - 1 e [-1/2, 1/2] 

are well defined and smooth with bounded derivatives on U{R) (as compo- 
sition of smooth functions with bounded derivatives). Moreover, we have 
ip — z{'ip)ri G W: as W is invariant under the multiplication by complex num- 
bers, the function 

n(^) := e-*"^'^^ - (1 + ri'ip))v = e-'^^^\i; - z{tP)rj) 

is in W, smooth for ip G U{R), and satisfies if) = x(a(V'))^(V');^('0))- 

To prove (4.5) let if)* E: T he the element of T realizing the minimum in 
the right-hand side (which exists by compactness of P). As ^* G F we have 
u{tp*) = 0. As the fonction tp i— )• u(V') is uniformly Lipschitz in U{R), we have 

||n(^)||^i <c||V'-V'*||^i=cdist(v^,r), 

which gives the result. D 

Let us now define the function u i— )■ r{u) from V^ to M by the implicit 
relation 

N{x{a,r{u),u)) = 1. 
By explicit calculation, we have 



(4.6) r(n) = -l + Vl-7V(n), 

from which we deduce that r{u) is well defined and smooth in a neighborhood 
of in H^, and moreover that ||r(ii)|| = 0(||n|| ^) if u is sufficiently small. 
Hence, (a, u) i— )■ xi<^j ^iu),u) is a local parametrization of 5 in a neighborhood 
of r C 5, where 

(4.7) S :={'tpeV\N{iP) = l}. 
Now let us define the function 

(4.8) n{u) = H{x{a,r{u),u)), 

which is well defined on W by gauge invariance of H. Moreover, this function 
is smooth in a neighborhood of 0. Then it can be shown (see [12]) that u = 
is a non degenerate minimum of Ti (u) : we have 



dniO) = 0, and VC/Giy, d^n{0){U,U) > c\\U\ 



H^ 
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Note that as % is smooth with locally bounded derivatives, the last coerciv- 
ity estimate extends to a neighborhood of uniformly: there exist positive 
constants c and p such that 

(4.9) VuG5(p), Vt/GW, d^n{u){U,U)>c\\U\\\^, 

where B{p) denotes the ball of radius p in W . In other words, the function % 
is strictly convex on B{p) and has a strict minimum at ti = 0. 

With these results at hand, let ^ G S, and assume that dist('0,r) is small 
enough so that we can write 

for some (a, u) £ T x W. Then for some constant C an sufficiently small u, 
we have 

dist(V',r) < 11^ - e'"r/||^, < Cir{u) + \\u\\^,) < C\\u\\^, . 

Now as n = is a minimum of the strictly convex function Ti on the ball B(p), 
we can write 

H{tlj)-H{ri) =n{u)-n{0) >l\\u\\\^ >cdist(V',r)2 

for some constants 7 and c > depending only on p. Then a Taylor expansion 
of 7i around n = shows that 

\n{u)-nm<c\\u\\\,, 

for some constant C depending on p and H but not on u € B{p). Hence using 
(4.5) we obtain the existence of constants c, C and Rq > such that for all 
ijj G S such that dist(V',r) < Rq, we have 

cdist(V',r)2 < \H{'^)-H{7])\ <Cdist('(/',r)2. 

The stability result (1.6) is then an easy consequence of this relation: As- 
sume that ipo £ S satifies dist(V'(0),r) < 6 < 5o where 60 < Rq, and let ip{t), 
t > be the solution of (1.1) starting at ■0(0) = V'o- Then by preservation of 
the energy H and norm A'^, we have 0(t) E S for all t > 0, and moreover as 
long as ip{t) is such that dist(V'(i),r) < Rq we can write 
(4.10) 
cdist(V'(t),r)2 < \H{iP{t)) - H{ri)\ = \H{%b{0)) - H{t])\ < Cdist(V'(0), r)^. 

Hence if 60 is small enough, this shows that for all t, dist(V'(t),r) < Rq and 
that (4.10) is in fact valid for all times t > 0. This implies (1.6) in the case 

N{i>) = 1. 



14 DARIO BAMBUSI, ERWAN FAOU & BENOIT GREBERT 



5. An abstract result 

In this section, we prove an abstract result for the existence and stabihty of 
discrete sohtons. We first give conditions ensuring that a discrete Hamiltonian 
acting on a discrete subspace of H^ possesses a minimizing sohton. We then 
show how the existence of a discrete flow (almost) preserving the Hamiltonian 
and the L^ norm ensures the numerical orbital stability over long times. In 
the next sections, we will apply this result to the three levels of discretization 
described above. 

5.1. Approximate problems. — We consider a set of parameter S € M^ 
and a function e : S — )■ R^ . This function will measure the "distance" between 
the discrete and continuous problems. 

For all ^ € S, we consider a Hilbert space Vfj, equipped with a norm || • || . For 
a given number R, we denote by B^{R) the ball of radius i? in V^. Moreover, 
for a given A; > a function F : V^ — )■ C of class C'', and a given ip^ ^ V^, we 
set for all n = 0, . . . , A; 

l|d»F(Vv)ll = .up |d"n^v)(C/'.....C/')l 



and we set 



sup sup ||d"F(V'^)|| 

n=0,...,fc iP^eB^{R) ' 



Moreover, we say that F is gauge invariant if it satisfies, for all a G T and all 
V'/i £ ^) F{e^"'il>^) = F{ipf^). Similarly, we say that G : V^ x V^ ^ C is gauge 
invariant if for all <p^ and V'/i in V^, and all a G T, we have G(e*"c/3^, e*"^^) = 

We assume that the family (V^);ies satisfies the following assumptions: 

(i) For all ;U G S, there exist a linear embedding i^i : V^ ^ H^ and a 
projection vr^ : H^ i— t- V^ that are gauge invariant in the sense that 
for all a G T and ijj^ G V^, e^^i^V/i = V^*""*/^/^ ^'^^ ^°^ ^^ V* ^ ^) 
e*"7r^^ = 7r^e*°V- Morever, we assume that i^ and tt^ are real in the 
sense that i^V^^ = i^V' ^'^^ '^^J.'4'|J. = '^f^'^i and that they satisfy the relation 
TT/i o V ~ ^"-^l Vm • Finally, we assume that there exists a constant Rq > 1 
such that for all /x G S, and c/?^ G B^{Rq), 

I 2 2 I 2 

(ii) For all yU G S, there exists a gauge invariant real scalar product ( • , • )^ 
such that setting Nf^{ip^) = {ip^,'>p^)^, we have Nf^{ip^) < \\ip^\\ and 
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(iii) For all ^ G S, there exists a gauge invariant function H^ : V^ — )■ M which 
is a modified Hamiltonian in the sense that 

(iv) If rj is the continuous soliton (1.4) defined in the previous section, we 
have for all ^ € S 

(5-1) Wit^^tiV-vWui < e(/^)- 

Note that using (i), there exist constants c, C and eo such that for ip^ £ V^ 
and n G T, such that e(/i) < eo, we have 

(5-2) cWi^ipJ^^ < UJ^ < Cllv^/'^ll^i . 

In the rest of this Section, we will assume that the hypothesis (i)— (iv) are 
satisfied. 

5.2. Local coordinate system. — We will assume here that all the /i G S 
considered satisfy the relation e(/x) < eo for some constant eo to be precised 
along the text. In echo to (4.7) we define for all // G S 

and the tangent space to tt^t] (compare (4.3)): 

W^ = {Uf, G V^ I {u^, iT^r])f, = {u^, ivr^??)/, = 0}. 

Note that if_iW^ is not included in W. 

By a slight abuse of notation, we will write u^ G -6/^(7) the ball of radius 
7 in W^ (instead of -6/^(7) n W^) for 7 > 0. We also set for i? > (compare 
(4.2)) 

(5.3) {V-M G V^ I dist^(V';„ ^;.r) < 7}, 

where dist^ denotes the distance measured in the norm || • || and where 



We then define the discrete application X/^ (see (4.4)): 

TxRxWf^B {a,r, u^) ^ Xf,{a, r, u^) = e*"((l + r)7r^r/ + u^) G V^. 

Lemma 5.1. — There exist constants eo, tq, C and R such that for all fi G 
S with e(/i) < eo, the application Xfi is smooth and bounded with uniformly 
bounded derivatives (with respect to n) from T x [— ro,ro] x B^{R) to V, and 
such for all ip^ G U^{R), there exists {a,r,Uf^) G T x M x 1^ such that ip^ = 
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Xfi{oi,r,u^). Moreover, the application Xn^ is smooth with uniformly bounded 
derivatives on U^{R), and for all ijj^ £ Ufj_{R), we have 

(5.4) \\uM^^)\\m <CdisV(V'^,7r^r). 

Proof. — The proof is exactly the same as the one of Lemma 4.3 by replacmg 
( • , • ) by ( • , • )^, A^ by Nfj_ and rj by n^rj. The fact that the constants are 
uniform in // is a consequence of the direct construction made in the proof of 
this Lemma and of the hypothesis (i)-(iv). Note that we use the fact that 

(5.5) IM^f^^) - 1| < Ce{i^)^ 

for some constant C independent on /x, which is a consequence of (ii) and 
(5.1), provided e(^) < eo is smah enough to ensure that ||7r^7?|| < Ri (which 

is possible upon using (5.1) and (5.2)). D 

Note that using the gauge invariance of i^, we have for all {a,r,u^) G 

T X M X W^ 

VX/i(a, r, u^) - x(a, r, i^u^) = e*"(l + r){i^TT^r] - rf) 
and hence for all u^ G W^, and r G M, 

(5.6) \\i^XtJ^{a,r,u^) - x{a,r,if,Uf,)\\^^ < (1 + |r|)e(^). 

Following the formalism of the previous section, we define for all ;U G S the 
function n^ i— )• r^{u^) on W^ by the implicit relation 

Nf,{x,i{a,rf,{Uf,),u^)) = 1, 

so that (a, u^) is a local coordinate system close to a rescaling of vr^F. Using 
the definition of A'^ and X/^j we immediately obtain that 



r^^{u^) = -1 + 






With this explicit expression, and using again (ii) and (5.6) there exist con- 
stants po, C and eo such that for all /i G S with e(/i) < eo, r^ is C'^{B^{pq)), 
and 

(5.7) II^M-'^°Vllc2(B^(po))<^e(/x), 

where the function r is defined in (4.6). Now defining (compare (4.8)) 

T-(-f^{Uf,) := H^{xi^{a,r^{Ufj),Ui,)), 

the previous relations, together with (iii) and (5.6) imply that if /?o is suffi- 
ciently small, T-L^ is well defined on B^{pq)^ and moreover 

(5.8) ll^°V-^.llc.(B,(po))^^^('^)- 

for some constant C independent of p, and for all ^u G S such that e(^) < eo- 
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5.3. Existence of a discrete soliton. — In the previous section, we liave 
shown that the continuous function % can be approximated by a function 
T-L^ on bahs of fixed radius pQ in Vs. This is the key argument to prove the 
following result: 

Theorem 5.2. — Under the previous hypothesis, there exists eo such that for 
all fi £ T, with e{pi,) < eo, there exists a discrete soliton r]^ £ V^ that realizes 
the minimum of H^ under the constraint N^{tp^) = 1, and such that 

(5-9) \K-^^ri\\y^<e{p). 

Moreover, there exist constants C , 5q and 70 such that for all fi £ Ti with 
e(/i) < eo, and all 5 < 5^, 

(5.10) dist(vVM>n' < C{\H^{il^^)-H^{7]^)\+e{f,) + 5), 
for all il)^ such that dist(i^'(/'/i)r) < 70 and \Nf^{ip^) — 1\ < S. 

Proof. — Let us take eo and po as in the previous section. Recall that as rj is 
a minimizer of the continuous Hamiltonian H, and by definition of 7^, we have 
d7^(0) = 0. Using (5.8), we deduce that for all /i G S such that e(/i) < eo, 

(5.11) ||d?^^(0)||^<C7e(;u). 

Moreover, for all U G W^, and n^ G B^(po), we have using again (5.8) 

\d^n^{u^){u,u)-d^n{i^u^){i^u,i^u)\ < Ce{fi)\\u\\l. 

Using (4.9) and (5.2), this shows that Tif^ is uniformly strictly convex in 
B/iipo), i.e. satisfies 

Vn^G5^(po), yUeW^, d^'H^,{u^){U,U)>co\\uf. 

with a constant cq independent on fi such that e(/i) < eo small enough. 

As Ti^ is strictly convex on the closed ball i3^(/0o)) ^^^l reaches its minimum 
on B^{pq) at some point u* G B^{po) (see for instance [6]). We want to prove 
that the minimum is reached in the interior of the ball. So assume on the 
contrary that u* is such that \\u*\\ = po, then we have 

^mK) - ^m(o) = d-H^m ■ n; + h{u;) 

with /i(u*) > co||n* II . Hence, as |d?^^(0) • n*| < Ce(/i)||u* || (see (5.11)) we 
get 

^m(^m) - ^m(O) > coPo - Ce{p)po. 
This shows that for eo sufficienly small, "Hfiiu*) > 'H^(O) which is a contradic- 
tion. Hence the u* is in the open ball B^(pq) and thus 

dn^iu;) = 0. 
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Moreover, as Ti^ is uniformly convex on the ball S^(/9o), we have 

\\u*\\^<C\\dn^{u*)-d-H^m^<Ce{fi). 
for some constant C independent on fj,. Then setting 

(5.12) r]f, :=X/.(0,r^(u*),'U*) = (1 + r^(u*))7r^r/ + n* , 

we verify using (5.7) and (5.12) that we have ||7r^?7 — r]^\\ < Ce{^) for some 
constant C independent on ^. 

It remains to prove (5.10). Let "0/^ S V^ and a € T, we have 

\\ • I 1 cv \\ ^ 11'/ 1 cy ' II ill" II 

where we used (5.2). Hence we have for all ij)^ 

(5.13) dist(v'(/'/.,r) < C7dist^(^^,^^r) + C7e(/i) 
for some constant independent of /i. Similarly we prove that 

(5.14) dist^(V^^,7r^r) < C7dist(v^^,r) + Ce(/i), 

for some constant C independent on //. Now let V^ be a function such that 
dist(i^VAi,r) < 70, with 70 small enough. Assume first that Nfj_{^^) = 1. 
Using (5.14), il)^ belongs to a set U^{'^) with a constant 7 depending on 70 
and eo- If these parameters are sufficiently small, we can define an element u^ 
of Bij,{pq) and a G T such that if)^ = X/i(0) '''fj,{ufj.),u^) (recall that Nf^{ip^) = 1) 
with u^ satisfying (5.4). Hence we have 

where u* is the minimizer of Ti, associated with the discrete soliton t/^. This 
implies that there exists a constant C independent of n such that 

\\u,-u;\\l<C\H,{i;^)-H,{r^^)\. 

Then using that ||n* || < Ce{fi), that ||u^|| = dist^(V'^,7r^r) + 0(e(^)), and 

the inequalities (5.13) and (5.14) we obtain (5.10) in the case N^{'iIj^) = 1. 
Now if Nfj^{'ip^) 7^ 1 but \N{ip^) — 1\ < 5 with 5 sufficiently small, there exists 
a point Vf^ such that ||f^|| < 6 and N{il)fj, — Vfj) = 1. We can then apply the 
previous estimate to ip^ — v^, and we use the uniform bounds on the derivative 
H^ to conclude. The approximation tp^ ~ v^ gives rise to the terms Cd in 
(5.10). D 
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5.4. Discrete orbital stability. — In the previous paragraph, we have 
shown that the conditions (i)— (iv) are sufficient to ensure the existence of a 
modified soliton for the modified energy //^, and that this soliton is sufficiently 
close to the exact soliton rj to allow the control of the distance between T and 
ipfj, via the distance between the Hamiltonian of H^{^n) and H^{r}^), see (5.10). 
As a consequence we obtain the following stability result 

Theorem 5.3. — Assume that the hypothesis (i)— (iv) are satisfied, and as- 
sume m.oreover that for all Rq and all fj, ^ Ti there exist /3{n) > and an 
application <1>^ : B^{Ro) — ^ V^ such that 

and 

(5.15) VV'^ E B^{Ro), \H^{^^{^l^^)) - H^{^^)\ < /3(/x). 

Then there exist 6o > and a constant C such that for all positive S < 5o and 
all fi £ T, such that e{fi) < eo and jp^ satisfying dist{if^ip^,T) < 6 then the 
sequence (V'^)n>o defined by 

vn>o, r^+^ = ^^ir^) 

satisfies 

Vn > 0, dist(vV';i, r) < C{S + e{^i)) 
as long as n/3{fj,) < e{fi) + 5. 

Proof. — Using the hypothesis on tp^ and (5.2), there exists Rq depending 
only on 6o such that Vu ^ B^{Rq/2) uniformly in // and there exists i^ € F 
such that Ili/iV'u ~ ^11 rri — ^' Thus using the gauge invariance of H, we have 
{Hii^tp^) — H{ri)\ < C6. Then with hypothesis (iii) and (5.9), we get 

|/f^(V^0)-F^(7?^)|<C(5 + 6(^)). 

On the other hand, using (5.15), we have for all n > 

n-l 

< C{S + e{n))+np{fi)<{C + l){6 + e{fi)) 
as long as nj5{iji) < e{iJ.) + 6 and HV'^H < Rq. Using the fact that Nfj,{ip^) = 
iV^(^O) = 1 + 0{6) and (5.10), we get 

(5.16) d[st{i^i^;,r)<C{6 + e{fi)) 

as long as llV'/ill < -Ro and for some constant C independent of /j, and n. 
Then by a bootstrap argument, there exists 6o and eo sufficiently small such 
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that, for < 5 < 6q and < e < eo, (5.16) ensures that this is the case for 
nj3{a) < e{fi) + 5. This proves the result. D 

6. Applications 

We now prove the three Theorems presented in Section 2. We only need to 
verify the hypothesis (i)-(iv) and to precise the constants e{fi) and /3(^). 

6.1. Discrete Schrodinger equation. — Consider the DNLS equation 
(2.2) for a given positive number h > 0. In the previous formalism, we set 
S = {/i E M"*"}, and the natural modified Hamiltonian and L^ norm are given 
by (2.3). We also define the real scalar product 

For all ^ G S, the embedding ih is defined by (2.5), and the projection vr/i by 
the application 

for some ip € H^. Defining the semi norm 

on V/j, we have by explicit calculation that 

(6.1) |V|^ = 14^1^1 

where \fp\ tti denotes the semi norm in H^. This fact allows to prove (i) and 
(ii) with the function e : S — )• M defined by e(/i) = h. This has already been 
proved in [2, Lemma 4.2]. Similarly, (iii) has been proved in [2, Proposition 
4.1] with e{h) = h. 

Finally, by classical arguments on finite elements approximation, there ex- 
ists an universal constant C such that for any function ip G H^ 

(6.2) \\7^h^P-^|;\\^,<Ch\m^,. 

This proves (iv) upon using (5.2). 

Let us define ^ji{ip) the fiow associated with the Hamiltonian Hh. Using 
standard estimates, one shows that it is well defined for sufficiently small t, 
say < i < to, uniformly in h. Theorem 2.1 is then a consequence of Theorem 
5.3 with 13(h) = and <1>^ = $^ with t G (0,to)- Remark that, in particular, 
since $"* = ($1)"" remains localized around the curve F of solitons for all n 
and for all t G {0,to), the flow $^("0) is deflned globally. 
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6.2. Dirichlet cut-off. — Recall that in comparison with the previous case, 
the space Vh^x defined in (2.7) is a finite dimensional space included in Vh- 
We have seen that the modified energy and norm Hh^K and Nh^x, and the 
embedding i^^K are defined by restriction to Vh^K- To define the projection 

T^h^K, we set 

r v(i/i) if IjI < K 

[ if \j\ > K. 

With these definitions, it is clear that the hypothesis (i)-(iii) are satisfied with 
S = {(/i, K) G M"*", xN} and with a priori e(/i) = h for /i = (/i, K). However, 
the estimate (5.1) is no longer true with the space cut-off without changing 
the definition of e{^). 

To have an estimate of \\ih,K''^h,K'n~v\\ iJi ) we only need to estimate \\T^h,K'n~ 
TT/jryll which is equal to 

\3\>K \j\>K 

By definition of ry, there exist constants Ci and v such that for all x G M, 
|r/(x)| < Cie~^\^K Substituting this estimate in the previous one, we get 

9 p — '^^jh 

\W,KV-^hn\\l < 2c!h Y, —^ + ^ E ^'^""^ 

\j\>K \j\>K 

n>K 

1 

< —K ex.p(-i'Kh) 

for some constant 7, and provided h < ho sufficiently small. 
This shows that (iv) is valid with the function 

(6.3) e(/x) = /i + -^exp(-zyis:/i), n = (h,K)£E 

With these notations. Theorem 2.2 is a consequence of Theorem 5.3 with 
/3(/.) = 0. 

6.3. Time splitting method. — Let us now consider the case where (2.6) 
is discretized in time by a splitting method of the form $^ o $p as described 
in Section 2. The space discretization being the same as in the previous Sec- 
tions, the hypothesis (i)-(iii) will be automatically fulfilled with the function 
e defined in (6.3). In particular, we can check directly that the norm Nh^x is 
preserved by splitting schemes. However, splitting methods do not preserve 
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the energy Hh,K for given h and K: more precisely, taking H^ = Hh^K in 
(5.15) only yields an error of order /3(/i) = hr. 
In this section, we set 

S := {{h,K,T) e M+ X N X M+}. 

For n = {h, K, t) G S, we set V^ = Vh,K, ifi. = ih,K = ih, and vr^ = tth^k- 
In the next section we will prove 

Theorem 6.1. — Let Rq > and M ^ N be fixed. There exist tq and ho 
such that for all // = (h, K,t) £ T, satisfying t < tq, h < ho and 

(6.4) (2M + 3)^<-. 

then there exist a constant C , depending only on Rq and M , and a smooth 
gauge invariant polynomial function H^ = Hhx,T defined on V^ such that 

(6.5) WH^-HoiJ^^^^^^^^^^Kc"^ 

and 

(6.6) ||<l>^ o cD^(V') - <^^hM)L < Ct^^' for all i; e V^ with \\i;\\ < Ro- 

With this result, the final statement of Theorem 2.4 is a consequence of 
Theorem 5.3 applied with 

1 r 

e(/i) = h + —K exjp(-uKh) + - 

n^ h 

and /3(/x) = r^^+^. The proof of Theorem 6.1 occupies the rest of this paper, 
and is a variant of the theory developed in [3, 11]. Here we summarize it 
and repeat the proofs with some details in order to have a quite self contained 
presentation. 



7. Construction of the modified energy 

7.1. Formal part. — We start by recalling the algorithm of construction of 
the modified energy Hfj^ introduced in the previous section. As a variant of the 
theory developed in [11], we work here at the level of the vector fields instead 
of Hamiltonian functions. Recall that at the continuous level, we identified 
the space H^(R;C) ~ if^(M;M)^ through the identification i/j = -ys{q + ip)- 
This identification obviously transfers to the space V of symmetric functions, 
and to the discretized space Vh^x via the identification 

(7.1) Vj = -j^iaj + iPj), J = -K, ...,K. 
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Hence we can endow Vh,K with the Hamiltonian structure induced by the 
symplectic form Yli=-K^Pj ^ ^^j- ^^ ^^^ following we make the constant 
identification between ■0 = (V'jO^-a- and {q,p) = {qj,Pj)f=^K given by (7.1). 
For a given real functional H[ip) = H{q,p), we associate the Hamiltonian 
vector fied Xh by 

ope oqi ) ^^_^ 

Note that this formula makes sense, because all the Hamiltonian functions 
-ff (V') that we consider are real valued. 

In this setting, A and P denote the vector fields associated respectively to 
the real Hamiltonian functions 

RA{i^) = h f; ^i^l^, and Hp{^) = -\ y: l^d^ 

e=-K e=-K 

which can obviously be expressed in terms of (qj,pj). Note that A and P 
depend on h, but we omit this dependence in the notation. We look for a 
formal vector field, namely a formal power series 

(7.3) Z(e):=J]Z,5", 

n>0 

where each Z„ is a Hamiltonian vector field on Vh^x, such that 

(7.4) V|e|<r, <P^p o <p\^ = <^l^^^, Ao := tA . 

Here <1>^ denotes the Hamiltonian flow on Vh^K associated with the vector field 

X at time t. 

Notice that, in particular, at order zero (7.4) implies 

(7.5) Zo := Ao = tA. 

Ideally, the approximate Hamiltonian we are looking for would be Hh,K,T '■= 
-Hz(r) (see (7.19)) but the formal series defining Z is not convergent and we 
will have to truncate the sum in (7.3). 

It is well known that it is convenient to look at the equality (7.4) in a dual way, 
namely to ask that the following equality is fulfilled for any smooth function 
w : Vh,K -^ C: 

(7.6) w{<^'po'^\^)=w{<^'zie))- 

The key ingredient of the construction is given by the formal formula 

(7.7) Vt, e*^^u; = 'u;o$*^ , 
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where Lx is the Lie operator associated with X. In our Hamiltonian case if 
X := {Xq,Xp)^^_j^ is a vector field (according to the decomposition (7.2)), 
we have in real coordinated {qj,Pj), 

and the exponential is defined in a formal way by 



1 

A;>0 






In this formalism (7.6) reads 

Deriving with respect to e one gets (by working on the power series) 

(7.8) e^^oe=^^Lp'u; = e^2(-')LQ(^)W, 
where 

(7.9) Q(£):=^^_^ad|(,)Z'(e) with adz X := [Z,X] , 

where [•,•] denote the Lie bracket of two vector fields. Finally (7.8) leads 
to the equation Q{e) = P from which we are going to construct Z{£). The 
construction goes as follows: first one remarks that the r.h.s. of (7.9) has 
the formal aspect of an operator applied to Z'{e), so the idea is first of all 
to invert such an operator. We remark that the power series defining the 
wanted operator is '^k>o^''/(^ + -'^)' ~ (^^ ~ 1)/^) so that one would expect 
its inverse to be x/{e^ — 1) = Ylk>o^''i-^k/kl), where B^ are the so called 
Bernoulli numbers and the power series is convergent provided |x| < 27r. So 
one is tempted to rewrite (5(e) = P in the form 

(7.10) VH<r, Z'(£) = J^|^ad|(,)P. 

k>0 

Plugging an Ansatz expansion Z{£) = J2i>o ^^'^i i^^to this equation, we get, 
for n > 0, the recursive equations 
(7.11) 



k\ 

fe>0 l.\^ |-4=n 



(n + l)Z.+i = j;^4'^\ with 4"):= Y. adz,^--.adz,^R 



Remark 1.1. — The analysis made to obtain this recursive equation is for- 
mal. To obtain our main result, we will verify that some of the series we 
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manipulate are in fact convergent series, while the others will be truncated in 
order to get nieaningfull expressions. 

Remark 7.2. — Assume that P is a polynomial of degree ro (in our case 
ro = 3), and that Zi is a collection of vector fields satisfying the previous 
relation, then for all n, Zn is a polynomial of degree (n — l)(ro — 1) + tq. 

Remark 7. 3. — If the vector fields P and Aq are Hamiltonian then the same 
is true for the vector fields Zn- This is an immediate consequence of the 
fact that all the construction involves only Lie Brackets, which are operations 
preserving the Hamiltonian nature of the vector fields. 

7.2. Analytic estimates. — We first introduce a suitable norm for mea- 
suring the size of the polynomials. In echo with the notations of the previous 
sections, we consider in the following a fixed /i = {h,K,T) € S. Recall that 
the space F^ = Vh,K does not depend on r, as well as the norm || • || . If 
X is a vector field on V^ which is a homogeneous polynomial of degree s 
we can associate to it a symmetric multilinear form Ar(^i, . . . , ipsi) such that 
X{iP) = X{i^,...,iP). We put 

ll-'^ll,, := sup ||X(V'i,...,V'si)||„ • 
'^ IIV'.I! =1 '^ 

M 

i=l,...,si 

We then extend this norm to general polynomial vector field X by defining its 
norm as the sum of the norms of the homogeneous components. 

Definition 7.4- — We denote by Vg the space of the polynomials of degree 
less than s, which furthermore have a finite norm || • || . 

Remark 7.5. — With this definition, we note that the norm ||P|| is uni- 
formly bounded with respect to ^. 

Lemma 7.6. — Let si > 1 and S2 > 1, and let X £ Vsi and Y G Vs2- Then 

[X,Y] CiVs^+s2-i, and 

(7.12) \\[X,Y]\\^<{s, + S2)\\X\\^\\Y\\^. 

Proof. — We give the proof in the case of homogeneous polynomials, the 
general case immediately follows. Denote again by X and Y the symmetric 
multilinear forms associated to X and y, then one has 

[x,Y]{^) = siX(y(v^),v...,V) - s2Y{x{^ij),^...,i^), 

from which the result immediately follows. D 
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Lemma 7.7. — For h < —^, the operator Aq = —rAh satisfies 



(7.13) Po||^<3^. 

Proof. — Let us first note that if {uj)^^_j^ is in V^, we have 

(7.14) Mi = 2h f; \^^-'p\\ h x; KP<(^+i)(/^ E K- 

j=~K j=-K j=-K 

Note that ^o = —T^f^ is homogeneous of degree one. Moreover, we can write 

(yloV'j^ = r 5 = -(o£ - a£_i), 

where ai = {ipe+i — ipi)/h. Using the discrete Sobolev inequaUty (7.14) and 
the Minkowski inequality, we get that 

K 



j=-K 



We conclude by remarking that 

K 



/^EKf)^' ' 



j=-K 

We deduce that 

\\Aoi;\\^<^V8 + 2h^\\ijJ , 

which shows the result. D 

Remark 7.8. — Lemmas 7.6 and 7.7 can be rephrased in a form suitable for 
the following by saying that, for X £ Vs, one has that the operator 

adx : Vsi -^ T's+si-i 
is bounded and its norm (induced by the norm || • || and for fixed s and si) 
fulfills 

(7.15) l|adX|| <(s + si)||X|| . 
\ / II II ^ — \ J-/ II II ^ 

In particular, using the previous result we have for a given si > 2 

T 

(7.16) ad^o : 'Psi -^ 'Psi and ||adAoll^ ^ 3(si + 1)-^ . 

Proposition 7.9. — Let M be an integer satisfying 

T 27r 



(7.17) (2M + 3)^< 



NUMERICAL SOLITONS FOR DNLS 27 

Then, for all n < M , Zn is well defined and Zn G Vr^ with rn = "^n + 1, and 
the norm of Zn is uniformly bounded with respect to /i. 

Proof. — We prove the proposition by induction. We set Zq = Aq. Assume 
that Zi S Vri for i < n < M — 1 are constructed. Let us prove that (7.11) 
defines a term Zn+i G ^r„+i- Rewrite (7.11) by incorporating the terms 
containing Zq = Aq and by substituting the estimate of the single terms to 
the ad terms. The advantage of doing that is that the product of the estimates 
is commutative, while the multiplication of the ad operators is not. We get 
first 

k 
II /(("•) II ^ Y^ II J ll'^"* ^' Y^ II J II II J II II Dll 

\\Al'\\ < > ad^n Tr TT—T > adz,, ... adz,, " 

i=l ^ ' ii + ...+ii=n 

"■ I. ■ y 

<y lladAoir ' ., \, , y (2r-M)"||%|| --WZe-W \\P\\ , 

j=l ^ ^ ii + ...+ii=n 

where we used that, if z > n and ij > then ii + ... + ii > n and the fact 
that, since by hypothesis the involved polynomials have degrees smaller then 
rM, one has || adz^ || < 2rjvf \\Zi\\ fj, foic i < n. 

Remarking that the result of the above sum with respect to £i , • • • ,£i does 
not depend on k, using (7.16) with si = ru, and noticing that ||P|| is uni- 
formly bounded with respect to fi, we get 

117 II ^ 1 ^^Bk^frurY'' kl ^ 
\\^-+^K ^ n + 1 ^ k\ ^ [ fi^ J (fc-i)!i!^" 

A;>0 i=0 ^ ^ ^ ^ ' 

~ n + 1 

for some constant C„ independent of /x. This shows that the series defin- 
ing Zn+i is convergent, that Zn+i G "Pr^+i and that ||Z„,+i|| is finite and 
uniformly bounded with respect to /i. D 

7.3. Proof of Theorem 6.1. — First remark that in our case all the vector 
fields are Hamiltonian. Explicitely, by Poincare Lemma, the Hamiltonian 
function of a Hamiltonian vector field X is given by 

(7.18) HxW:= I suj{X{s^),^)ds , 

Jo 

where u is the symplectic form. In particular, this formula shows that the 

Hamiltonian function of a smooth polynomial vector field is also a smooth 
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polynomial function. For e < r, let us define 

M 
j=0 

By construction Z^^'^'{e) satisfies (7.10) up to order e^ included from which 
we deduce that it satisfies (7.4) up to order e (see [11] Theorem 4.2 for 
details). Therefore defining for // = {h,K,T), 

1 ^ 

(7.19) H, := -Hz(M)(,) = Y^r^-'Hz^ , 

j=0 

estimate (6.6) holds true with a constant independent of /i. 

It remains to compare the two Hamiltonians Hf^ = Hj^x^t and Hy^x in the 
C^ norm on the ball centered at the origin and of arbitrary radius i?o in V^. 
Let us define 



t'^B. 



k\ 



' adi, P . 



and recall that Zq = Aq = tA, and that by construction 

(7.20) ^i = E 

fc>0 

Now we have 

M 

i=2 
But using (7.18) and the fact that Zj is of degree rj, we get for all ip G i?^j(i?o)) 

M 

(7.21) |i/«(^) - H,m < J2r'''\\Zj\\^Rl^^' < Ct , 

i=2 

for some constant C independent of h, K and t < tq sufficiently small. To 
estimate iJ^ — Hh^x, we notice using (7.20), 



(7.22) 



^rrW - XHh,K = Zl- P = T 



Y^i^^d^ 



k>0 



{k + iy. 



Ao 



ad^o P- 



But in view of (6.4), ^ < vr, and thus the operator ^^^q (k+i)\ ^"^Ao 



IS 



bounded on Vs, uniformly with respect to fi. Therefore for ip £ B^{Ro), we 
have 

i4i)(v)-^/.,if(v)i<c|F,d^^p(V')i. 
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for some constant C independent on /i = (/i, K, r). Now we calculate explicitly 
that the Hamiltonian associated with adyio ^ i^ given by 

K 



-K 
K 



T 



^9((V'^+i + V<?-i -2V'^)l^£pV'^ 



But we have 

K 



i=-K 
K 



=-K 
K K-1 



e=-K i=^K-i 

K-l 



i=-K 

using the boundary conditions V'x+i = i'-K-i = 0. Taking the imaginary 
part, we obtain 

K-l 

HadAgPW = ^ ^ ^{{A+i-i^e){\i^efi'e-\'4'e+i\^i^e+i))- 

But we have 

|9((V'<?+i - 'ipe){\tpi\'^ipe - \ipi+i\'^ip£+i))\ < 5\ipe+i - •i/'£p(|V'£p + IV'^+il^)- 
Then we use that 

to obtain 

|i:^ad^„pWI<2-||V||' 



T 



and therefore for ip £ i3^(i?o) 

(7.23) \H^^\i;)-Hh,Km<C 
Combining (7.21) and (7.23) we get, for all ip G B^{Rq), 

(7.24) \H^{ilj)-Hh,Km<C- . 
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Furthermore, since both functionals are analytic in ip and the above estimate 
is uniform in ^ G i3^(i?o)) we have similar estimates for the first and the 
second derivative of ^ i— )■ H^{tp) — Hh^Kiip)- C 
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